Non-Perturbative Functional Renormalization Group for Random Field Models and 
Related Disordered Systems. II: Results for the Random Field O(N) Model 
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We study the critical behavior and phase diagram of the d-dimensional random field O(N) model 
by means of the nonperturbative functional renormalization group approach presented in the pre- 
ceding paper. We show that the dimensional reduction predictions, obtained from conventional 
perturbation theory, break down below a critical dimension (1dr{N) and we provide a description 
of criticality, ferromagnetic ordering and quasi-long range order in the whole (N, d) plane. Below 
dDn(N), our formalism gives access to both the typical behavior of the system, controlled by zero- 
temperature fixed points with a nonanalytic dimensionless effective action, and to the physics of 
rare low-energy excitations ("droplets"), described at nonzero temperature by the rounding of the 
nonanalyticity in a thermal boundary layer. 

PACS numbers: ll.10.Hi, 75.40.Cx 



I. INTRODUCTION 

The random field model describes one of the simplest 
disordered systems in which classical ./V-component vari- 
ables (spins in magnetic language) with O(N) symmetric 
interactions are linearly coupled to a random (magnetic) 
field. Yet, more than thirty years after the first studies of 
the modelji^ its long-distance behavior (criticality and 
ordering) still largely appears as a puzzle. 

The two main questions raised about the (equilibrium) 
properties concern the nature and the characteristics of 
the phases and of the phase transitions. The first one 
is about the so-called "dimensional reduction" property 
which relates the critical behavior of the random field 
O(N) model (RFO(iV)M) in dimension d to that of 
the pure O(N) model in dimension d — 2. This prop- 
erty, predicted to all orders by conventional perturba- 
tion theory^ii and derived as a consequence of a hidden 
supersymmetry^ is known to break down in low enough 
dimensions. -i2 The second question concerns the exis- 
tence of a phase with quasi-long range order (QLRO) 
(i.e., a phase characterized by no magnetization and a 
power-law decrease of the correlation fuctions) in the 
models with a continuous symmetry (N > 1) below 
d = 4, their lower critical dimension for long-range ferro- 
magnetism. More specifically, the presence of QLRO in 
the 3-dimensional KFXYM (N = 2) is of relevance to the 
"Bragg glass" phase discussed in the context of vortices 
in disordered type-// superconductors i 8 : 9 : 10 ' 11 : 12 

In the preceding article,™ denoted as paper I in the 
following, we have developed a nonperturbative functional 
renormalization group (NP-FRG) formalism to study the 
long-distance physics of random field models and related 
disordered systems. In the present article, we put the 
formalism to use to address the issues mentioned above 
concerning the behavior of the RFO(iV)M. 

We start in section II by briefly recalling the main def- 
initions, notations and results of the NP-FRG approach 
to the RFO(iV)M which have been presented in paper I. 



Next, in section III, we discuss the mechanism by 
which dimensional reduction breaks down, namely the 
appearance of a strong enough nonanalytic behavior in 
the field dependence of the dimensionless effective aver- 
age action. We first recall the analysis of the perturbative 
functional RG at one loop near d = 4. We then present 
the scenario for the failure of dimensional reduction in 
the RFIM within our NP-FRG approach and extend our 
considerations to the whole (N, d) plane. 

The numerical results obtained from our minimal tru- 
cation of the NP-FRG for the RFO(iV)M are presented 
in section IV. This allows us to provide a unified de- 
scription of criticality, fcrromagnetism, and QLRO in the 
whole N — d diagram. We show that two nontrivial crit- 
ical lines characterize the long-distance behavior of the 
RFO(iV)M (on top of the upper, d uc = 6, and lower, 
dic(N = I) = 2 and di c (N > 1) = 4, critical dimensions 
for the paramagnetic to ferromagnetic transition): a line 
dDn(N) separating a region of the (N, d) plane in which 
the critical exponents are given by the dimensional re- 
duction predictions (d > dr)n(N)) from a region where 
dimensional reduction is fully broken (d < doniN)); and 
a line di c (N) characterizing the lower critical dimension 
for quasi-long range order (for 1 < N < N c — 2.83...). 
Finally, we discuss the accuracy and reliability of the 
present truncation. 

In section V, we address the physical meaning of the 
nonanalyticity (in the effective action) which is associ- 
ated with dimensional reduction failure. We also exam- 
ine the role of temperature and the connection with the 
phenomenological droplet approach. In particular, we 
discuss the "activated dynamic scaling" behavior that 
characterizes the critical slowing down of relaxation in 
the RFIM. 

We finally conclude by considering the relation of the 
present to other pictures of the behavior of random field 
models and providing some perspectives. 

Short accounts of the p resent work have already been 
published in Refs. 
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II. NP-FRG APPROACH FOR THE RFO(AT)M 

The NP-FRG approach developed in paper I combines 
three main ingredients: 

1. a version of Wilson's continous RG in which one 
follows the evolution of the effective average action 
Tfc with a (momentum) scale fc, from the bare ac- 
tion at the microscopic scale (fc = A) to the full 
effective action at macroscopic scale (fc = 0); the 
evolution of T^, which is the generating functional 
of the 1-particle irreducible vertices at scale k, is 
governed by an exact RG flow equation. 16 

2. A replica formalism in which the permutational 
symmetry among replicas is explicitly broken by 
the introduction of linear sources acting indepen- 
dently on each replica; using an expansion in the 
number of unconstrained (or "free") replica sums 
gives access to a description of the probability dis- 
tribution of the renormalized disorder through its 
cumulants. 

3. A nonperturbative approximation scheme for the 
effective average action that relies on truncating 
both it "derivative expansion" (expansion in the 
number of spatial derivatives of the fundamental 
fields) and the "expansion in number of free replica 
sums" (or, equivalently, the cumulant expansion). 

For the RFO(AT)M, the minimal truncation of 
which already contains the key features for a nonpertur- 
bative study of the long-distance physics is the following: 

~ s n n 
Tfe [{0a}] = / \- 5>m,fc|^ a W| 2 + £ U k (<t> a ( X )) 
Jx ^ a=l o=l 
1 " 

- 5 E V*(0.(s),&(*)) 



a, 6=1 



(1) 



where as before 4> a , a = l,...,n, are the replica fields, 
Z mt k is a wave function renormalization parameter, Uk 
is the 1-replica potential which physically represents a 
coarse-grained Gibbs free energy and gives access to the 
thermodynamics of the system, and Vk is the 2-replica 
potential which is the second cumulant of the renormal- 
ized disorder evaluated for uniform fields. 

The flow equations for Uk(<j>i), Vk{4>\,4>2), and Z m ^ 
are obtained from the exact RG equation for the effective 
average action ; 13 ! 16 



d k T k [{<t> a }} 



Tr\d k R k (q 2 



,(2) 



Ri 



. (2) 

where the trace involves a sum over both replica indices 

(2) 

and A- vector components and T k is the tensor formed 
by the second functional derivatives of T k with respect 
to the fields 4>a{<l)- Rk(q 2 ) is the infrared cutoff which 



enforces the decoupling of the low- and high- momentum 
modes at the scale k. It is diagonal in A^-vector indices, 
and in the minimal truncation we have also^ chosen it di- 
agonal in replica indices, i.e., R k " ab {q 2 ) = Rk(q 2 )S a bS IJ , v . 

The initial condition is given by the bare replicated 
action of the RFO(iV)M, 

S[{4> a }}= I {^J2[\d ( f>a^)\ 2 +r\<f>a^)\ 2 + 

J<* < A1 =l 

a, 6=1 ' 

(3) 

where we have made explicit the dependence on a "bare" 
temperature T. 

One more step is needed to cast the NP-FRG 
flow equations in a form suitable for searching for 
the anticipated zero-temperature fixed points of the 
RFO(A r )M ) 17 i 18 namely, to introduce appropriate scaling 
dimensions. This requires to define a renormalized tem- 
perature Tk which is expected to flow to zero as k — > 0. 
Near a zero-temperature fixed point, one has the follow- 
ing scaling dimensions: 



Tk ~ k , Zn 



fc" 



it - kH d - i+f >\ 



(4) 



with 9 and rj related through 8 = 2 + r\ — r], as well as 

U k ~ k d ~ e , V k ~ k d - 26 , (5) 

so that the second cumulant of the renormalized random 
field, 



Af (<j> u <p 2 )=d<j>1d<p%V k (<j>i,cl> 2 ), 



(6) 



scales as fc ' 2 '' v \ (We recall that the superscripts with 
greek letters denote the components of the JV-vcctor 
fields.) 

The dimensionless counterparts of U k ,Vk, Afe,0 are 
denoted by lower-case letters, Uk,Vk,Sk,<p. (For the 
RFO(A r )M, a convenient parametrization of the 1- and 
2-replica functions makes use of the variables p — |y| 2 /2 
and z = ipi ■ ip 2 / '\fApxp~2.) The resulting flow equations 
in scaled form have been given in section IV of paper I. 

We conclude this brief recapitulation of the results of 
paper I by recalling that the minimal nonperturbative 
truncation described above reduces to the 1-loop pertur- 
bative results near the upper critical dimension d — 6 and 
when A^ — > oo, and, most importantly, that it reproduces 
the 1-loop perturbative FRG equations near the lower 
critical dimension for ferromagnetism in the 0(N > 1) 
model, d = 4. 



III. BREAKDOWN OF DIMENSIONAL 
REDUCTION 

As stressed in the introduction, dimensional reduc- 
tion for the random field model must break down in low 
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enough dimension. Within the NP-FRG, we find that 
the mechanism by which this occurs is the appearance 
along the RG flow of a nonanalyticity in the field depen- 
dence of the dimensionless effective average action; more 
precisely, the appearance of a cusp in the second cumu- 
lant of the renormalizcd random field as one makes the 
two field arguments approach each other. The theory is 
renormalizable, albeit with the unusual feature that the 
renormalized effective action is nonanalytic at the fixed 
point. Such a mechanism has previously been found in 
the random manifold model within the perturbative FRG 
approach^* 2 ^ 2 ^ 

In this section, we discuss the appearance of a non- 
analytic behavior within our minimal truncation scheme. 
We first recall the results obtained near d = 4 for the 
RF V(N > l)M^i 24 i 25 i 26 since this limit is more easily 
accessible to an analytic treatment and already provides 
the scenario for the general case. 



A. RFO(A)M at one loop near D = 4 

Our starting point is the set of equations derived at 
first order in e = d — 4 and at zero temperature for 
the running exponents rjk, fj k , an d for the dimensionless 
renormalized second cumulant of the disorder, more pre- 
cisely for R k {z) = v k (p mt k,p m ,k,z)/(2pm,k) 2 where p m ,k 
corresponds to the minimum of the 1-replica potential 
(i.e., is akin to a dimensionless order parameter at the 
running scale k) and goes as 1/e at the relevant fixed 
points: see Eqs. (100) and (101) of paper I. (Note that 
Rk(z) should not be confused with the regulator Rfc(q 2 ) 
that appears for instance in Eq. |(3J).) For studying the 
fixed points and their stability it is convenient to intro- 
duce Rk(z) = (ivi/e)Rk(z) and to rescale the RG "time" 
as et — > t. The flow equation then reads 



d t R k (z) = R k (z) - 2(N - 2)R' k {\)R k {z) - l(N 



1) 



R' k (z) - 2zR' k (l)) R' k {z) - -(1 - z 2 ) [ - R' k {zf 



2(i4(l) - zR' k {z))R' k \z) + (1 - z 2 )R' k \zf 



where t = ln(k/A), and one also has 

r, k =eR' k (l), f} k =e[(N- l)R k (l) - 1], 



(7) 



(8) 



It is instructive to start by analyzing the flow equa- 
tions for the first derivatives R' k (z = 1) and R k (z = 1), 
assuming that Rk(z) is at least twice continuously diffcr- 
entiable around z = 1. these equations read: 



-d t R' k (l) = -R' k (l) + (N-2)R' k (l) 2 



-«(1) = (-l+6i4U)K'(l) 



+ (N + 7)RUl) 2 + R' k (l) 2 



(9) 



(10) 



If Rk(z) is analytic around z = 1, the flow equations 
for the higher derivatives evaluated in z = 1 can be de- 
rived as well. As noted by Fisher,— the expression for 
the pth derivative only involves derivatives of lower or 
equal order. This structure allows an iterative solution 
of the fixed-point equations obtained by setting the left- 
hand sides to zero, provided of course that R*(z) has the 
required analytic property. 

Beside the stable fixed point R*(l) = 0, there is one 
nontrivial fixed point associated with Eq. (|9j): 



K(l) = l/(N-2), 



(11) 



with a positive eigenvalue Ai — e. This fixed point leads 
to the dimensional-reduction value of the critical expo- 
nents, i.e., i] — fj — e/(N ~ 2), v — 1/Ai = 1/e. On the 
other hand, Eq. (jTUJ) has nontrivial fixed-point solutions 
only when TV > 18. These solutions are: 



{N-8) + ^/{N~2)(N-18) 
2(N-2)(N + 7) 



(12) 



which is unstable with an eigenvalue A2 
and 



(JV-18) 
(JV-2) 



(N — 8) — y/(N -2)(N -18) 



2(N -2)(N + 7) 



(13) 



which is stable with A2 = 



jw^ 6 - For N < i 8 > no 

fixed-point solutions exist for Eq. pop . One instead finds 
that there is a (finite) range of initial conditions R'Al) 
for which the RG flow for R k (l) leads to a divergence at 
a finite scale k, irrespective of the initial value R'^(l). 

The solution to the absence of a nontrivial, twice diffcr- 
entiable fixed-point function R*(z) when N < 18 is that 
the proper fixed point controlling the critical behavior is 
nonanalytic around z — 1, with R'^(z) having a cusp, i.e., 
a term proportional to \fl — z when z — > 1. Numerical 
solutions showing this cuspy behavior have been given by 
Fcldman 25 for N = 3, 4, 5 and by us for general values of 
N < 183 

We have shown in detail in Ref. [26| that the value 
Ndr = 18 separates a region in which R'^z) at the criti- 
cal, i.e., once unstable, fixed point has a cusp (N < Ndr) 
from a region (iV > Ndr) where R'„(z) has only a weaker 
nonanalyticity, a "subcusp" in (1 — 2) a W with a(N) 
a noninteger larger than 3/2. The occurence of a cusp 
changes the values of 77 and fj from the dimensional re- 
duction prediction, t\dr = Vdr = e /(-^~ 2). Indeed, the 
flow equation for R' k (l) is modified according to: 



-d t R' k (l) = -R' k (l) + (N-2)R> k (l) 2 + 

lim {2(1 - z)R' k \z) [2(1 - z)R' k "{z) - ZR'^z) 



(N + l)R' k \z)~ 2(1 -z)R'l'(z) (R' k (z 



#*(!))}. 
(14) 
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where the whole term lim z ^i(...) is nonzero when a cusp 
is present in R' k {z). As a result, the once unstable fixed- 
point solution for R' t (l) is no longer equal to 1/(TV — 2) 
and it follows from Eqs. ijHJ) that dimensional reduction 
is broken^ 

On the other hand, the weaker nonanalyticity occuring 
for TV > 18 does not alter the flow equation for R' k (l), 
which is still given by Eq. ([9]) , and dimensional reduction 
still applies; in particular, rj — fj — tjdr- This drastic 
change of behavior at TV — 18 is illustrated in Figure Q] 
where t^dr/v and fjoR/fj are plotted as a function of TV. 




FIG. 1: Breakdown of dimensional reduction near d = 4: 
Vdr/ti (lower curve) and fjoii/fj (upper curve) versus TV for 
the RFO(TV)M with TV > 3 at first order in e = d - 4. The 
dimensional reduction value of the exponents is t)dr = Vdr = 
e/(N-2). 

To get some insight into the order of the nonanalyticity, 
one may analyze the hierarchy of flow equations for the 
successive derivatives of Rk(z) evaluated in z = L 24 i 26 As 
explained above, a fixed point with a well defined second 
derivative and an associated negative eigenvalue (A 2 < 0) 
can be found for TV > 18 (see Eq. ([T5|)^ 

Let us assume that the first p derivatives of Rk (z) are 
well defined in z — 1. Contrary to the flow equations for 
R\(l) and Rf k \l) (see Eqs. (|9TT0)) ). those for the higher 
derivatives are linear, namely, 



-d t R { k p) (l) =A p (R! k (l),R!ai)) i?i p) (l) 
+ TF p (i4(l)X'(l), 



,i?i p - 1} (l)), 



(15) 



where A p and T v are known functions easily derived from 
Eq. ||7J). If R' k {l) and R k (l) are chosen equal to their 
fixed-point values given in Eqs. (|llll3p . one finds that 



Ap * ~ TV - 



2p 2 -(TV-l)p+(TV-2) + 



p(TV - 5 + 6p) 
2(TV + 7) 



(TV - 8 - v/(TV- 2)(TV- 18)) 



(16) 



For a given TV, there exists an integer value p$(N) such 
that A p » < forp < pj(TV) and A p * > forp > p i (N) + l. 



The RG flow for the (pj (TV) + l)th derivative therefore di- 
verges when t — > — oo whereas all lower-order derivatives 
reach finite fixed-point values. As a consequence, the 
fixed-point function R'^(z) must have a nonanalyticity of 
the form (1 - z) a W with p B (TV) - 1 < a(TV) < pj|(TV). 
Refining the reasoning, — one finds that a(TV) is given by 
the solution of A Q (jv)+i* = 0, where A^^+i* is given 
by Eq. (fTrj|) with p replaced by the noninteger a (TV) + 1. 
The result is shown in Figure [2j the order of the non- 
analyticity increases with TV when TV > 18, starting from 
3/2 when TV -> 18+, and it goes as TV/2 + 0(1) at large 
TV. 



a 




16 18 20 22 24 26 28 30 

TV 

FIG. 2: Exponent a(TV) characterizing the order of the non- 
analyticity in the cumulant of the renormalized random field, 
R'(z), near z = 1 for the RFO(TV)M in d = 4 + e. A cusp 
leading to breakdown of dimensional reduction corresponds 
to a(TV) = 1/2 and is obtained for TV < 18. There is a dis- 
continuity at TV = 18 with a(TV = 18+) = 3/2. 



Finally, we stress the different ways in which "cusp" 
and "subcusps" appear along the RG flow. As seen 
above, subcusps occur only at infinite RG time (i.e., at 
the fixed point). On the contrary, due to the nonlinear 
nature of the beta function for R k (l) (see Eq. (jlOp ). a 
cusp appears at a finite time, which one may define as a 
"Larkin scale" by analogy with the behavior of disordered 
elastic systemsi 22 ' 23 i 27 ' 28 



B. Analytic versus nonanalytic behavior in the 
RFIM 



Consider now the Ising version of the random field 
model at zero temperature. We have seen in section IV-B 
of paper I that "anomalous" contributions to the com- 
bination of running exponents 2r/k — r)k may appear if 
the dimcnsionless renormalized cumulant 5k(^>i,f2) be- 
comes nonanalytic as the two field arguments approach 
each other, ipi — > ip\. This can be more conveniently 
studied by changing the variables to x = (tpi + </?2)/2 
and y = (ipi — ip2)/2; the dependence of 8 k on x is an- 
ticipated as being completely regular and that on y as 
potentially anomalous. (Recall that due to the Zi and 
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permutational symmetries, 5 k is even in x and y sepa- spect to ip\ and tf2 and switching to the new variables x 
rately.) The flow equation for Sk(x,y) is obtained from and y. One finds that 
that for v k {ip±, ^2) given in paper I by deriving with re- 



d t 5 k {x,y) = (2rj- rj)6(x,y) + hd - 4 + rj){xd x +yd y )S(x,y) - l { 2 d \w+)S 0+ + l {d) (w-)5 - 

- 2l[ d l(w + ,w.)S(x,y)]s^(x,y) - l[ d l(w + ,w^ 01 Hx,y) 2 + [(4^+)^+ - if (w-)S' _) 
+ 2{l {d l{w + ,w^)w' + - lW(w + ,w-)v/_)8{x,y) - 2(l (d) (w+)w' + 5 Q+ - lf(w^w'_5 ^]s^(x,y) 
+ [i d \w + )6 0+ -lW(w-)6 ^ n \x,y) + 2l^ 

+ [d d) (w + )S' 0+ + li d \w-)6' _) 2(li d \w + )w' + 6 0+ + i d \w-)w'_S Q -) - 2(li d l(w + ,w.)w' + 
+ (w + , w _ )w'_ )S(x, y)] 6™ (x, y) + \ [l {d) (w + )S Q+ + l {d) («,_)*<>- + 2l (d } (w + , W -)S(x, y)] <5 (20) (x, y) } , 

I 



(17) 



where we have dropped the subscript k in the right- 
hand side and introduced the short-hand notation w± — 
u k( x ± V)> s o± = $k,o(x ± y) with 5kfi(x) = 5 k (x,y = 
0); for a function of a single argument a prime de- 
notes a derivative whereas for functions of two argu- 
ments partial derivatives are indicated as superscripts, 
e.g., 6W(x,y) = d x 5{x,y), S^(x,y) = d y 5(x,y), etc. 

Finally, In (w) and ii?,n 2 (wi,W2) are the "dimension- 
less threshold functions" defined from the infrared cutoff 
function R k {q 2 ) = Zkq 2 r(q 2 /k 2 ) (see paper I and Ref. 

We now follow a reasoning similar to that developed 
for the RFO(iV)M near d = 4. Assume that Sk{x,y) is 
continuously differentiable with respect to y around y = 
up to some order 2p with p > 1. Then, introducing the 
notation 5 k , q {x) = d!*5 k (x, y)\ y =o and using the property 
that all derivatives of odd order vanish in y = due to 
the inversion symmetry, one may express 5 k (x,y) in the 
vicinity of y = as 



and for the derivatives, 



5 k {x,y)=j2i^5k M {x)+o{y 2 n. 



q=Q 



(25)! 



(18) 



By inserting this expression in the RG flow equation for 
5k(x,y), Eq. (fT7|) . one derives the following flow equa- 
tions for the function evaluated in y = 0, 



d t Sk,o( x ) = ( 2r lk - Vk)Sk,a(x) + 4 + Vk)x8' kfi (x)- 

2v d {l (d \uUx))S ki0 (x) 2 u^(xf-4 d) (u'l(x)) u k "(x)x 

SkAxWkA*) + 4%*(s))[§**,o(*) a + s kt oWU*)]}, 

(19) 



d t S k ,2{x) = -L 2 [u k ,S kt0 ]Sk,2(x) + Zv d l 2 d) {u k \x))5k,2{x) 2 

- 2v d g 2 [u k ,8 k fi], 

(20) 



and for p > 2, 
dtS k ,2 P (x) = -L 2p {u k ,8 kl o, 5 k ,2]Sk,2 P (x) 

- 2v d g 2 p[u'l, {5 k ,2q} 



(21) 



g<p-lJ' 



where L 2 [u k ,5k,Q] and L 2p [u k , 6 ki o, 5k, 2] are linear oper- 
ators whose expressions are given in Appendix A. The 
C?2p's are functionals of u' k (x), 5 k ,o(x), and of the deriva- 
tives 5 k ,2 q (x) with q < p — 1. Their expressions are not 
worth displaying here. 

The above equations are complemented by the flow 
equation for u k (x), or its derivative u' k (x), obtained from 
the results of section IV-B in paper I, 

d t u' k {x) = -(2 - r)k)u' k (x) + i(d - 4 + f} k )xu k (x)+ 

2v d {l[ d \u' k \x))5' k ^x) - l < i\ul{x))u'l!{x)5 k>0 {x)}, 

(22) 

and the expression for the running anomalous 
dimension, 13 



T]k = y|2m^(«j(3; m i),^(i m i.))iif(i mi t)' 
m 2 d 2( u k\xm,k),u k (x mt k))u k (x m ,k)\, 



(23) 



where x m ,k denotes the nontrivial configuration that 
minimizes the 1-replica potential, and which therefore 
satisfies u' k (x m k ) — 0; the ml?^ n2 (wi, w 2 Ys are addi- 
tional "dimensionless threshold functions" (see paper I 
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and Ref. [l6|]). Similarly, one also has an expression for 
2% ^ fjk which is derived from Eqs. (Til?]) and (|22|) and 
the constraint 5k,o{xm,k) = 1 (one recovers the equation 
of paper I, but now without the anomalous terms). 

One first notices that the RG equations for u' k (x), 
8k,o{x), and rjk form a closed set. No aditional input is 
required from the derivatives 5k,2p{x) with p > 1, which 
means that the RG flow for the 1-replica potential and 
for the 2-replica potential (or the second cumulant) eval- 
uated for equal field arguments is closed without further 
knowledge of the full field dependence of the 2-replica 
potential for distinct replicas. This property is a direct 
consequence of the assumptation that the behavior of the 
second cumulant is sufficiently regular when the two field 
arguments become equal, more precisely, that dy8k(x,y) 
is finite when y — > 0. 

Before discussing the consequences of this property, it 
is worth mentioning that it results from the structure of 
the exact RG equations and not from the specific approxi- 
mation chosen here. More generally indeed, the exact RG 
flows for the 1-replica component of the effective average 
action and for the cumulants of the renormalized random 
field (see sections II-C,D of paper I) evaluated for equal 
field arguments decouple from the full functional depen- 
dence of the cumulants when the latter is regular enough 
in the limit of equal arguments. This point will be fur- 
ther developed and clarified in a forthcoming publication 
centered on the superfield formalism^ 

As in the previously discussed case of the RFO(7V)M 
near d = 4, one expects that the fixed point obtained 
without reference to distinct replicas, and associated with 
a regular enough behavior of the cumulants in the limit of 
equal arguments, corresponds to dimensional reduction. 
To prove this, one needs to show that it is equivalent, in 
the 1-replica sector at least, to the corresponding fixed 
point of the pure system in two dimensions less. This 
is indeed illustrated near the upper critical dimension 
d = 6: it is easy to show (see also section IV-D of paper 
I) that, at first order in e = 6 - d, Eqs. (TT51 1221 |35|) 
give back the result of the pure Ising model at first order 
in e = 4 — d. The difficulty in going beyond this step 
is that the present truncation does not necessarily pre- 
serve the underlying supersymmetry of the model. For 
instance, the second order in e = 6 — d of Eqs. (THil [221 
I23p breaks the dimensional reduction property; however, 
this is clearly an artefact of the truncation and of the 
choice of regulator. One can check this by improving the 
treatment so that the exact 2-loop results are recovered 
near d = 6: the calculation becomes extremely tedious 
in the present formalism and a nonperturbative closure 
becomes hardly tractable numerically; but it is nonethe- 
less found by a direct analysis near d = 6 that if the 
2- and 3-rcplica cumulants are regular enough in their 
field dependence so that the flow equations evaluated for 
equal replica fields decouples as in Eqs. (fH)l 1221 123"j) . the 
corresponding fixed point at second order in e leads to 
dimensional reduction. 



Awaiting for a proper resolution of the problem via 
the superfield formalism^ we will associate with dimen- 
sional reduction the fixed point corresponding to Eqs. 
(flUl [221 123"]) . fixed point that can be continuously followed 
as a function of dimension d and, as will be discussed 
below, as a function of the number of components N. 
Breaking of dimensional reduction therefore implies the 
occurence of a strong enough nonanalyticity in the field 
dependence of the renormalized cumulants of the disor- 
der. "Strong enough" here means that it is sufficient to 
couple the flow of the components of the effective aver- 
age action evaluated for equal fields to the full functional 
dependence involving distinct replica fields. 

From the flow equation for 6k(<pi,<P2) = 6k(x,y), Eq. 
(1171) . it is clear that the only way to avoid dimensional 
reduction is therefore the existence of a linear cusp in the 
fixed-point function 5<,(x,y), i.e., with 

6,(x,y) = 8^ (x) + \y\5^ a (x) + O(y 2 ), (24) 

in the vicinity of y = 0. This leads to the appearance 
of an "anomalous" contribution to the expression of the 
beta function for 6k,o(x): 

P So (x) = f3 So \re g (x) - V d lf «'(x))5*, (z) 2 , (25) 

where — (3$ \ reg (x) is given by the right-hand side of Eq. 
(Til?]) . The above equation gives back the expression of 
2rjk — fjk derived in section IV-B of paper I with the tem- 
perature set to zero. 

Before closing this discussion of the RFIM, we would 
like to emphasize a few additional points which parallel 
the comments made in the previous subsection. First, 
the appearance of a cusp in 8k(x,y) is associated with 
the divergence of the second derivative 8k,2{x). Due to 
the nonlinear character of the flow equation for dk,2 (x), 
Eq. (j20|) . this divergence, if present, is expected to first 
occur at a finite scale which, as before, we generically call 
the "Larkin scale" . For a running scale k larger than the 
Larkin scale k^, the effective average action is analytic 
and it develops a cusp (in 8k(x,y)) for k less than kh- 
From the flow equation for 5k(x,y), one can see that a 
linear cusp is stable under RG flow, in that it does not 
lead to stronger "supercusps" . 

Secondly, the appearance of a "subcusp" in Sk(x,y) is 
signaled by the divergence of a higher-order derivative 
8k,2p(x), with p > 2 being related to the order of the 
nonanalyticity (which is strictly less than 2p and strictly 
more than 2(p— 1)). The flow equation for 5k,2 P {x) being 
linear, a subcusp can only appear at infinite RG time, 
i.e., at the fixed point. Following the reasoning devel- 
oped in Appendix A, we conclude that the order of the 
nonanalyticity characterizing the fixed point increases as 
1/(6 — d) 2 when d approaches 6 from below. This is the 
counterpart of the situation found near d — 4 where the 
order of the nonanalyticity increases as N/2 when the 
number of components N gets large. In both cases the 
critical fixed points with fully analytic (dimensionless) 
effective action found above d = 6 and when N — > oo 
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are approached in d or N by fixed points with diverging 
orders of the nonanalyticity, i.e., with weaker and weaker 
subcusps. We recall however that such subcusps are not 
sufficient to break dimensional reduction. 



C. Extension to the whole (N,d) plane 

The preceding developments on the connection be- 
tween breakdown of dimensional reduction and nonan- 
alyticity of the effective average action can be extended 
to the whole (N, d) plane. The nonanalyticity now oc- 
curs in the renormalized cumulants of the disorder as two 
arguments, i.e., two replica fields, approach each other, 
V2 — * <Pi- If the nonanalyticity is weak enough, namely, 
if it is weaker than a linear cusp in the second cumu- 
lant of the renormalized random field, the RG flows for 
Uk (p) A-,r (p) ,5k,L(p), Vk (see section IV-C of paper I) de- 
couple from those involving distinct replica fields. The 
associated fixed point can be continuously followed in 
the (N, d) plane and, near d = 6, near d = 4 for N > 18 
and when N — > oo, it leads to dimensional reduction. 
For reasons explained above, a direct proof that it corre- 
sponds to dimensional reduction away from the pertur- 
bative regime is hampered by the truncation used here, 
but we rely on the continuity argument within the (JV, d) 
plane to nonetheless identify it with an approximation of 
the dimensional reduction fixed point. 

Breaking of dimensional reduction is thus associated 
with the appearance of a sufficiently strong nonanalyt- 
icity in the 2-replica potential Vk(pi, pi, z). Analysis of 
the flow equation for Vk{pi, P2, z) (see paper I) shows that 
this corresponds to the following behavior as <Pi—*(pi: 

Vk(pi,P2,z) ^ v k , reg (pi, P2, s 2 ) + \s\ 3 v k , a (p,r 2 , \s\), (26) 

with s 2 = \ipi-ip 2 \ 2 /{8y/pTp2), p = {pi+p 2 )/2, r = {pi- 
p 2 )/(A\s\), and \ipi-ip 2 \ 2 = 2(pi + p 2 - 2^/p 1 p 2 z); v k , r eg 
and v kyll are analytic functions of their arguments in the 
vicinity of p 2 = p\ = p, z = 1 {i.e., s = 0), and r 2 < p 2 . 
The cusp in the second cumulant of the renormalized 
random field, S k v (pi,p 2 , z) = d^d^Vkipi, p 2 , z), occurs 
in \s\ or, equivalently, in \tpi — tp%\; it is marked by the 
divergence of the second derivative of v k with respect to 
s 2 when s = (which also implies p 2 — pi). Note how- 
ever that on top of \<£>i — ip 2 \, there is now an additional 
variable, denoted r above, that characterizes the way tp 2 
approaches tp\. 

To conclude this section, it should be stressed that 
the consistency of the present scenario and the actual 
occurence of a cusp in a region of the (N, d) diagram 
must be verified by a numerical resolution of the NP- 
FRG equations. This is what we address now. 



IV. A UNIFIED DESCRIPTION OF 
CRITIC ALITY, FERROMAGNETISM AND QLRO 

A. RG flow equations and their numerical 
resolution 

The RG flow equations for the RFO(A^)M in the mini- 
mal truncation of the NP-FRG discussed above are given 
in paper I. Focusing on the fixed points and their vicinity, 
we drop the subdominant terms involving the tempera- 
ture (it will be checked that the temperature exponent 6 
is indeed strictly positive). The structure of the resulting 
equations can be summarized as follows: 

dt.u' k (p) = -f3 u >[u' k ,6 k , T ,Sk,L;Vk,Vk}(p), (27) 
d t v k (pi,P2,z) = -f3 v [u' k ,v k ;r] k ,fik}{pi, p 2 , z), (28) 

Vk = lri(Pm,k, u'k{pm,k),h,T(Pm,k),h,L{pm,k)), (29) 

where (3 u i and f3 v are functionals and j ri is a func- 
tion; S k ,T(p) = (2p)~ 1 d z v k (p,p,z)\ z=1 and S kiL (p) = 
2pd pi d P2 v k (pi, p 2 , z = l)\ pi=P2=p are the transverse and 
longitudinal components of the second cumulant of the 
renormalized random field evaluated for equal field ar- 
guments, and p m , k is the configuration that minimizes 
the 1-replica potential (u k (p m , k ) = 0). The running 
exponent fj k is derived from the flow of the constraint 
&k,T{pm,k) — 1 that follows from the definition of the 
renormalized temperature^ 

Eqs. (|27H29|) form a set of coupled partial differen- 
tial equations involving in particular a function of three 
variables. Studying the whole (N,d) plane by numeri- 
cally solving these equations remains a very difficult and 
computationally intensive task. To facilitate the study, 
we have used in addition an expansion in powers of the 
fields. However, some caution must be exerted in order 
to retain enough of the functional character for allow- 
ing a description of possible cusp or nonanalytic depen- 
dence. We have therefore considered an expansion of 
pi,p 2 around the configuration p rn . k while keeping the 
complete dependence on the variable \(p\ — tp 2 \ 2 which 
we anticipate to be the key variable for describing the 
cusp (see above). More specifically, we have chosen the 
following approximation: 

«*(P)= (¥) (P-P^k) 2 , (30) 



v k {pi,p 2 ,z) = v kfi0 (s 2 ) + v ktW (s 2 )(pi + p 2 - 2p m , k )+ 

7^fc,2o(s 2 )(/n + P2 - 2p m ,k) 2 + 7^fc,02(s 2 )(pi - p 2 ) 2 , 

(31) 

wh ere s 2 = \(p 1 - <p 2 \ 2 / {8p m ^ k ) = (pi + p 2 - 
2yJpip 2 z)/(Ap m ^ k ) (when p x = p 2 = p m . k , s 2 varies 
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then between and 1); by construction, v' koo (0) — 
Sk,T(Pm,k) = 1- It is easily checked that this trunca- 
tion still reproduces the perturbative results near d = 6, 
d = 4, and when N — > oo (compare with section V of 
paper I). 

Inserting the above expressions into the flow equa- 
tions, Eqs. H27H29p . provides a set of cou- 
pled partial differential equations for 4 functions, 
Vkm(s 2 ),Vkfii{s 2 ),v k ,2o{s 2 ),Vkm(s 2 ), and 3 running pa- 
rameters, Pm,k-> ^k,Vk (phis, when convenient, fjk), whose 
resolution now represents a more tractable numerical 
problem. 

We close this subsection by outlining the numerical 
methods used for solving the partial differential equa- 
tions. For each couple (N, d) and for a choice of the 
infrared cutoff function (in most calculations, we have 
taken the "optimized" regulator 30 described in paper I 
that leads to explicit analytic expressions for all the di- 
mensionless threshold functions appearing in the beta 
functions), we follow the evolution under RG flow of the 
various functions and parameters for given initial condi- 
tions. For the functions, a finite difference mesh is used 
for the variable s, so that standard algorithms are suf- 
ficient to solve the evolution with the RG "time" t (or 
equivalently, the scale k). When necessary, most notably 
for checking the robustness of a nonanalytic cusp-like be- 
havior around s = 0, we vary the mesh spacing. To reach 
the critical, once unstable, fixed point, we fine tune the 
initial condition for p m k which represents the unstable 
direction. We consider that a fixed point is attained when 
the sum of the absolute values of all beta functions is less 
than 10~ 6 . The whole procedure can be accelerated by 
following the fixed points by continuity (when possible) 
in the (N, d) plane through small finite changes of N 
and/or d. 

We now move on to the presentation of the main re- 
sults. 



B. Dimensional reduction and its breaking 

To study the "weakly nonanalytic" critical fixed point 
associated with dimensional reduction (see section III), 
one does not need the full dependence on s of the func- 
tions Vk,oo, Vk.oi, v k.2Q, Vfe,02) but only their value and that 
of their first derivative evaluated in s = 0. In addition, 
to check the stability of this fixed point to the appear- 
ance of a cusp in |s| in the dimcnsionless cumulant of the 
renormalized random field Sk (or, equivalently, a term in 
|s| 3 in Vk), we also follow the second derivative of the 
functions, evaluated in s = 0. 

We find that the dimensional reduction fixed point 
is stable versus cusp-like behavior in a whole region of 
the (N, d) plane. However, there is a critical dimension 
<1dr(N) depending on N (which one may as well de- 
scribe for fixed d, as a critical number of components 
NDR.{d)) at which the second derivative of Vk with re- 
spect to s 2 in s = first diverges along the RG flow at 



a finite "Larkin" scale. The difference in behavior above 
and below dp^N) is illustrated in Figure [3] For N = 3, 
we display the evolution with t of the second derivative 
of Vk with respect to s 2 evaluated in p\ = p2 = Pm.k 
and s — 0, i.e., up to a constant prefactor, v'^ 00 (0). 
The initial condition on p m ,k has been fine tuned so that 
the other running quantities (see above) reach the di- 
mensional reduction fixed point. For d = 5.5, vju OO (0) 
reaches a finite fixed-point value; conversely, for d — 5.0, 
it diverges at a finite Larkin time. (Note that the other 
second derivatives v'l 10 (0) , v'l 20 (0) , v% 02 (0) all diverge at 
the same Larkin scale.) The value of dDn(N) in this case 
in about 5.1. 



V 




1 2 3 4 5 6 

-t 

FIG. 3: Difference in behavior above and below cIdr{N) for 
N = 3: evolution with RG "time" \t\ of v'l 00 (0) which, up to 
a constant prefactor, is the second derivative of the dimen- 
sionless 2-replica potential Vk with respect to s 2 oc \ipi — (p2\ 2 
evaluated in pi — p2 = p m .k and s = 0. The initial condi- 
tion on Pm,k has been fine tuned so that the other running 
quantities reach the dimensional reduction fixed point. For 
d = 5.5, Vfc'oo(O) reaches a finite fixed-point value; conversely, 
for d = 5.0, it diverges at a finite "Larkin" time. 

When d — > 4 + , we numerically recover the value 
Ndr = 18, thereby confirming that the nonperturbative 
truncation actually leads back to the perturbative FRG 
result near d = 4. The curve NjjR(d) extends continu- 
ously down to N = 1, where we obtain duR(N) ~ 5. It 
separates the two regions denoted I and IV in Figure S] 

For d < rfrjjj(iV), the fixed point controlling the crit- 
ical behavior of the RFO(7V)M must now be studied 
by keeping the full dependence on s of the renormal- 
ized disorder cumulant. We do find in this region that 
a cusp occurs at a finite time that corresponds to the 
Larkin scale discussed above. In Figure [5] we illustrate 
the change of behavior related to the presence or ab- 
sence of cusp in the fixed-point function <5*T(Pm,fc, s 2 ) for 
N = 3. We plot 6*T(Pm,k, s 2 ) as a function of both s 2 
and d. Below some dimension 5.1, a cusp is clearly visible 
near s — 0. We have checked its robustness and that it 
is indeed a behavior in |s| by changing the spacing of the 
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d 

FIG. 4: Predicted phase diagram of the d-dimensional 
RFO(iV)M. In region III, there are no phase transitions and 
the system is always disordered (paramagnetic). In regions 
/ and IV, there is a second-order paramagnetic to ferro- 
magnetic transition and in region II, a second-order transi- 
tion between paramagnetic and QLRO phases. In region IV 
the nonanalyticity of the dimensionless effective action at the 
zero-temperature fixed point is weak enough to let the critical 
exponents take their dimensional reduction value, whereas a 
complete breakdown of dimensional reduction occurs in re- 
gions / and II. 




FIG. 5: Presence or absence of a cusp in the fixed-point 
function 5*T(pm,k, s 2 ) for N = 3: 2-dimensional plot of 
<5*T(Pm,fc 5 s 2 ) versus s 2 and d. Below some dimension around 
5, a cusp is clearly visible near s = 0. 



discretization of the s variable. The value of d for which 
the cusp first appears coincide with cLdr as determined 
previously from the divergence of d^ 2 Vk\ P = Pm k s =o] here, 
d DR {N = 3) ~ 5.1. 

Finally, we draw attention to the fact that the critical 
exponents evolve continuously upon crossing the critical 
line dDn{N). This has been clearly illustrated in Fig- 
ure [2] for the exponents 77 and fj near d = 4: one can see 
that fj gradually separates from 77 as TV moves down from 
18 and does not settle in general to the value fj = 2-q. 
This provides strong evidence in favor of a characteriza- 
tion of the critical scaling behavior by three i ndep endent 
exponents and not two as suggested in Refs. [31| . 



C. Criticality, ferromagnetism, and QLRO 

In the present NP-FRG formalism, the equilibrium 
phases of the system are characterized by the flow evo- 
lution of the dimensionless order parameter at scale k, 
p m ,k- Depending on the initial conditions, p m ,k is found 
to:' 

1. diverge when time goes to —00 (k — > 0) in such 
a way that the dimcnsionfull order parameter (the 
magnetization) tends to a finite value; this behav- 
ior can be associated with a fully stable (attractive) 
fixed point describing long-range ferromagnetic or- 
der, 

2. go to zero at a finite scale, which corresponds to 
the disordered paramagnetic phased 

3. reach a nontrivial fixed point value; the dimension- 
ful order parameter which goes as k d ~ 4+r, p m * is 
then zero when k = 0. 

In the latter case, if only one parameter (say, the initial 
value of p mt k) need be adjusted, the fixed point is once 
unstable and describes the critical point of the model 
and the associated scaling behavior. This is the situation 
found in regions I and IV of the (N, d) diagram shown in 
Figure[4]with, as discussed above, a qualitative difference 
with respect to dimensional reduction between the two 
regions. 

A different pattern is found below dimension 4 for mod- 
els with continuous symmetry (N > 1). Whatever the 
initial conditions, the RG flow no longer leads to a diver- 
gence of p m ,k, which is in line with the predicted absence 
of long-range order in the RFO(iV)M for d < 4. One 
obtains instead a nontrivial attractive fixed point char- 
acterized by a cusp in the second cumulant of the renor- 
malized random field. This fixed point comes on top of 
the once unstable (critical) fixed point also characterized 
by a cusp. It describes a whole low-disorder phase that 
is associated with a nontrivial scaling behavior. This 
phase has a vanishing order parameter but algebraically 
decaying correlation functions characterized by the two 
anomalous dimensions 77 and fj (see section IV-A of pa- 
per I). It therefore corresponds to a QLRO phase and it 
transforms into the disordered paramagnetic phase at a 
critical point, itself controlled by a "cuspy" fixed point. 

The QLRO phase only occurs below a critical value 
of the number of components, N c = 2.83..., which can 
also be directly computed from an analysis of the per- 
turabative FRG equations at and near d = 4. 15 i 26 ' 32 It is 
very similar to the QLRO phase found in elastic systems 
pinned by disorder ^iSJi and for the RFXYM (N = 2) 
it actually identifies to the latter phase when d — + 4~. 
(The numerical solution of the truncated NP-FRG equa- 
tions then reduces, as expected, to the 1-loop perturba- 
tive FRG result for which the equivalence with a random 
periodic elastic model is easily shown: see e.g. paper 
I.) However, contrary to the situation in disordered elas- 
tic systems, the QLRO phase only exists below a (bare) 
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critical disorder and the existence of two "cuspy" fixed 
points provides a mechanism for destroying QLRO below 
some lower critical dimension . 15 i 26 i 32 

The variation with N and d of the characteristics of 
the two "cuspy" fixed points is illustrated in Figure 6 
where we plot the value of the anomalous dimension r\ 
versus d for a series of values of N. For N > N c ~ 
2.83, only one fixed-point value emerges from the point 
(7/ = 0,d = 4); but for N < N c , one finds two values 
of 77 for each dimension, the upper one being associated 
with the critical fixed point and the lower one with the 
QLRO fixed point. One can see that the two fixed-point 
branches coalesce for a value di c (N) which consequently 
determines the lower critical dimension below which no 
phase transition is observed. 




FIG. 6: Anomalous dimension n associated with the single or 
the two zero-temperature "cuspy" fixed points plotted versus 
d for values of N ranging from 1.4 to 4 by steps of 0.2. For 
N > N c ~ 2.83, only one fixed-point value emerges from the 
point (77 = 0, d = 4); but for N < N c , one finds two values of 
77 for each dimension, the upper one being associated with the 
critical fixed point and the lower one with the QLRO fixed 
point. The two branches of fixed points coalesce for a value 
di c (N) shown by filled circles. 

The lower critical dimension is shown in Figure |4j 
where it separates region II with QLRO and region III 
with no phase transition. Note that the result is compat- 
ible with Feldman's prediction^ 3 , that there is no QLRO 
above N — 3. The lower critical curve, seen as Ni c (d) 
in place of di c (N), decreases as d decreases and is ex- 
pected to reach iV; c = 1 for d = 2, which corresponds 
to the lower critical dimension (for long-range ferromag- 
netism) in the Ising version. (In our present approximate 
description, we find Ni c (d = 2) ~ 1.15 instead of the ex- 
act result Nic = 1: see the discussion below.) This be- 
havior is in fact reminiscent of what occurs in the pure 
O(N) model.— Although never acknowledged, the (N, d) 
phase diagram of the latter, derived from the solution of 
phenomenological RG equations^! and from known exact 
results, is indeed very similar to that of the RFO(iV)M 
displayed in Figure S) the critical value N c is now equal 
to 2 and occurs in d = 2, which is the lower critical di- 
mension for ferromagnetism; a QLRO phase exists in a 



region between (N = l,d = 1) and (N = 2,d = 2), 
which is the counterpart of region II described above, 
but is of course not physically relevant for real systems. 
This QLRO is different from that associated with the 
low-temperature phase of the XY model in d = 2, be- 
low the Kosterlitz-Thouless transition ) 35 i 36 as the latter 
corresponds to a line of fixed points while the former is 
controlled by a single fixed point. 

The situation found in the RFO(iV > 1)M below d = 4 
raises two questions. First, is a QLRO phase observ- 
able in an experimentally realizable random field sys- 
tem ? The question is especially acute because it has 
been suggested that the 3 — d RFAFM displays such a 
phase ^ 10 ' 37 called "Bragg glass" in relation with pinned 
vortex lattices in disordered type-// superconductors, 
the second question is more academic and concerns the 
status of the singular point (N c = 2.83..., d — 4). Fol- 
lowing the analogy with its counterpart, (N c = 2, d = 2), 
in the pure O(N) model (see above), one may wonder 
whether the random field model with (N c = 2.83... in 
d = 4 also gives rise to a Kosterlitz-Thouless transition. 
Elsewhere ; 15 ' 26 we have answered this question through 
an analysis of the perturbative FRG equations near d = 4 
at 2 loops, finding that there is no Kosterlitz-Thouless- 
like transition. A practical consequence of this result is 
that, as indeed reproduced by the numerical solution of 
the truncated NP-FRG equations, the critical line Ni c (d) 
drops abruptly with an infinite slope as one moves away 
from d = 4 (see Figure [4J . 

Going back now to the first question raised above, 
we conclude from looking at Figure 2] that the N = 2 
RFA"Y"M in d = 3 is very distinctly in the region where 
no QLRO phase occurs. Therefore, we find that there 
is no Bragg glass in this model. This is confirmed by a 
direct study of our equations in JV = 2, d = 3J£ As a 
cautionary note here, we would like to point out that the 
absence of a Bragg glass phase in the 3 — d RFAYM does 
not necessarily implies that no such phase exists in dis- 
ordered type-// superconductors. The relation between 
vortices in the latter systems and the RFAVM is derived 
via an "elastic glass model" ^ and it does not guaran- 
tee that phase transitions associated with the presence of 
massive modes are identical in the two systems. 



D. Robustness of the results 

The results presented above rely on a nonperturbative, 
but approximate, RG description. It is therefore desir- 
able to have some estimate of its degree of accuracy. In 
addressing this point, we make use of two main proper- 
ties of the present NP-FRG formalism, which we have 
emphasized in several occasions in this article and in pa- 
per I: 

1) the existence of a systematic truncation scheme that 
allows one to control and improve the results by going to 
higher orders of the truncation; this has been tested with 
success on the pure O(N) mode l 16 ! 38 ' 39 , and we plan to 
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do in the near future for the RFO(A)M (the calculations 
being, however, much more demanding). 

2) a unified description of the whole (N, d) plane for 
the d-dimcnsional RFO(A)M; one can therefore check 
the consistency of the approximate nonperturbative re- 
sults by comparing with known exact or perturbative re- 
sults in the appropriate regions of the (N, d) plane. This 
is what we address in the following. 

The numerical resolution of the truncated NP-FRG 
equations confirm the property already stressed above 
and in paper I, that one recovers the 1-loop perturba- 
tive predictions near d = 6, when N — > oo, and, more 
significantly, near d — 4. A study of the 2-loop pertur- 
bative FRG equation s 15 ! 26 ' 32 near d = 4 unambiguously 
supports the scenario found here concerning the pres- 
ence of two nontrivial critical lines, one associated with 
the breaking of dimensional reduction, N DR (d), starting 
downward from N — 18 when d — > 4 + , and one giv- 
ing the lower critical dimension for QLRO below d = 4, 
Nic(d), arriving in N c = 2.83... with an infinite slope 
when d — > 4 _ . There are some quantitative differences 
(see Ref. [26j]), but the overall picture near d = 4 is well 
captured by the present truncation. 

To estimate the error made in locating the critical line 
NDR.(d), we have also considered a somewhat cruder ap- 
proximation for the 2-replica potential Vk, which is a 
plain expansion in powers of the fields (still around the 
minimum) including all terms up to Lp A : 

Vk = 2«l,fe(VPlP22 " Pm,k) + V 2 ,k(Pl + 92 - 2pm,kf 
+ V 3tk (Pl - Pi) 2 + V4,fc(V/°l/°22 - Pm.k) 2 
+ V 5 . k {y/piP2Z - Pm.k)(Pl + P2 ~ 2Pm,k), 

(32) 

where Vi k = 1 by construction. The results (obtained 
from monitoring the divergence of d? Vi _ ip , 2 Vk\ v , 1 = 1 p 2 , see 
above) are very similar to those obtained with the more 
involved truncation: the critical line Ndr (d) starts again 
from TV = 18 near d = 4 and extends down to don, — 5 
when N = 1. The value of doniN = 1) is found in the 
window 4.9 — 5.1, depending on the approximation and 
the choice of cutoff function.— 

Finally, one may also compare the lower critical value 
Ni c (d) obtained when d = 2 to the expected exact re- 
sult, Nic(d = 2) = 1. As stated above, we find Ni c (d = 
2) ~ 1.15, which provides an estimate of the error. (The 
lower critical dimension of the RFIM is a difficult an un- 
favorable test for the present approach starting from a 
Ginzburg-Landau- Wilson bare action; as shown for the 
pure Ising model, the results can be improved by solving 
the full first order of the derivative expansion^) 

The above discussion indicates that there is certainly 
room for improvement of the quantitative predictions 
(and we have provided a formalism to do so), but it 
also gives strong confidence in the robustness of the 
present description of the long-distance physics of the 
RFO(iV)M. 



V. PHYSICS OF THE CUSP AND ROLE OF 
TEMPERATURE 

We devote this section to a discussion of the physics 
associated with the nonanalyticity found in the effective 
average action and of the role of temperature. To pro- 
ceed, we build upon the body of work already done in 
the context of disordered elastic system o 20 ' 21 ' 41 ' 42 ' 43 . For 
ease of notation, we focus on the RFIM in the follow- 
ing, but most results will apply mutatis mutandis to the 
RFO(A)M. 



A. Interpreting the cusp 

Breakdown of dimensional reduction has been associ- 
ated with the presence of multiple "metastable states" , a 
metastable state being generically taken as a field config- 
uration that minimizes some action or effective action. 
From the supersymmetric formalism of the RFIM^i 4 ^ 
one finds that a necessary condition for this breakdown 
is the existence of many minima of the bare action (as 
given by Eq. (1) of paper I). In the present approach 
on the other hand, the failure of dimensional reduction 
predictions originates in a strong enough nonanalyticity 
in the field dependence of the dimensionless effective av- 
erage action. To shed light on the connection between 
this nonanaliticity and a picture in terms of metastable 
states, we follow the line of reasoning developed for elas- 
tic systems pinned by a random potentiah 20 ' 21 ' 41 ' 42 ' 43 

To begin with, it is worth stressing the unusual char- 
acter of the RG analysis in the presence of a cusp. Gen- 
erally speaking, integration over fluctuations, e.g., ther- 
mal fluctuations in Statistical Physics, smooth away non- 
analyticities as well as the effect of possible metastable 
states so that, at long distance, the dimensionless effec- 
tive average action is a nonsingular function of the fields. 
The novelty in the RFIM case comes from the dominance 
of the (quenched) disorder fluctuations over the thermal 
ones and the associated property that the long-distance 
physics (criticality, ordering and quasi-ordering) is con- 
trolled by zero-temperature fixed points. Actually, this 
physics is describable by working directly at zero temper- 
ature at all scales: see above and paper I. As argued in 
the context of disordered elastic systems / 9 ' 20 ' 21 integra- 
tion over high-energy modes along the RG flow amounts 
at zero temperature to minimizing some coarse-grained 
action, and it is this minimization procedure that may 
lead to cusp-like behavior in the presence of multiple min- 
ima in the coarse-grained action. 

It is now instructive to go back to the interpretation of 
the 2-replica potential Vk{4>i, ^2) as the second cumulant 
of the renormalized disorder and of its second derivative 
Afc (</>!, 4>2) as the second cumulant of the renormalized 
random field, both being evaluated for uniform field con- 
figurations (see paper I). To be more precise, the 1-replica 
component of the effective average action r^ i^] is the 
Legendre transform of the first moment of the random 
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free energy functional at scale k, Wk[J; h], namely, 

r M [<f>] = -W k [J; h] + f J(x) ■ 4>(x), (33) 

where J(x) is a linear source conjugate to the field <j>(x) 
and the overbar denotes an average over quenched disor- 
der, with h denoting the bare random field. On the other 
hand, the 2-replica component is the second cumulant of 
W k [J; h] with J = Jfc [</>] , where J k [<ft] is the nonrandom 
source defined via the above Legendre transform, i.e., 
Jk[<j>](x) = 8Tk.i[4>]/ $4>{x). One therefore has 



r fe)2 [<k, H = 8W k [J k [<h]-,h]5Wk[Jk[<h]-,h], (34) 



with 5W k [J; h] = W k [J; h] - W k [J- h}. 

One can also define a renormalized random field at the 
running scale k, h k [(/)](x), as 

hMx)=--^SW k [J[4>]-,h]. (35) 

It has zero mean and its second cumulant is given by 



h[<h]{x)hk[<h](v) = T k% y l<f>i,fa}- (36) 

More details, including a discussion of higher-order cu- 
mulants, can be found in paper I. 

In the truncated NP-FRG considered here, the random 
free energy functional <5Wfc[J[</>]; h] is taken in a local ap- 
proximation which amounts to replacing it by a random 
potential V k (<p;x) with zero mean and second cumulant 



% [{<t>x;x)V k (<fe ;y)^6(x-y)V k ( ( f >1 ,h)- (37) 

Similarly, in this approximation, the renormalized ran- 
dom field defined above is given by h k (4>]x) = 
—d<j>V k [(4>i',x) with a second cumulant 

h k {4>i\x)h k {4> 2 ;y) ^ 6(x - y)& k (<t>i,fa)- (38) 

All the above considerations of course apply to the 
dimensionless quantities, v k (f\,ip2) and 5 k (ipi,tp2). In 
particular, one can introduce a dimensionless random po- 
tential iik with its second cumulant given by v k (tfx,f2) 
and an associated dimensionless random field with sec- 
ond cumulant given by 5 k (tpx,ip2)- In what follows, we 
rather discuss the dimensionless functions since, at the 
fixed point, the dimensionful quantity V k goes to zero 
whereas Afc diverges. 

Following Ref. [2l| . a cusp in <5*(<^i, if 2) as <f% — > <f\ can 
be interpreted as resulting from a cuspy random potential. 
Such a potential is sketched in Figure 7a. It displays a se- 
quence of minima separated by cusps located at random 
positions along the field axis. (The explicit spatial depen- 
dence of £>* (f, x) has been dropped as for instance could 
be obtained from a properly rescaled integration over 
space.) As seen in Figure 7b, the dimensionless random 
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FIG. 7: Sketch of the field dependence of the dimensionless 
renormalized random potential v*(ip) (a) and random field 
h*(tf) = — dtpV*(<p) (b) associated with a cusp in 6*((pi,tp2) 
as f2 <pi. The cusps separating the minima in (a) corre- 
spond to dicontinuities or "shocks" in (b); they are located at 
random positions along the field axis. 

field then shows discontinuities at those random loca- 
tions., discontinuitites that can be tnough of as "shocks" 
through an analogy with the Burgers equatio m 21 i 45 

To see how the rugged cuspy landscape of Figure 7a 
gives rise to a cusp in the second cumulant of the random 

field 6*((pi,(p2), consider the quantity (h^fi) — h*((f 2 )) 2 
as f 2 - * fi- Switching again to the variables x — (ip\ + 
y>2)/2 and y = (<f\ — <p-i)l% one has from Eq. (f38|) : 

(K(fi) - K(tf 2 )) 2 ^ 2(5* (x,0) - 6*(x,y)) , 
~-25*, a (x)\y\+0(y 2 ), 

where we have used the cusp-like behavior as y — * (see 
section III-B) and where 5*, a (x) should be negative. 

On the other hand, with a renormalized random field 
as pictured in Figure 7b, (h*((fi) — h^(f2)) 2 is a 0(y 2 ) 
except when a discontinuity (shock) is present between 
the two fields f± and if 2- As a consequence, 

— z rx+\v\ r 

(h*(<pi) - K(if 2 )) 2 ^ / dx d dj d p(xd,Jdhd 

J X -\y\ J 

+ 0(y 2 ), 

(40) 

where Xd is the shock location and 7^ the ampli tude of t he 
associated discontinuity. Assuming that J drfdP(x,jd) r Yd 
is nonzero, one indeed recovers Eq. ([39]) with S*. a (x) = 

- J d ld p(x^d)l 2 d - 

The above discussion therefore points to a picture in 
which, when dimensional reduction is broken, the fixed 
point controlling the critical behavior of the RFIM is de- 
scribed (for uniform field configurations) by a dimension- 
less random potential with multiple minima separated by 
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cusps. These minima and the cuspy barriers separating 
them arrive along the RG flow, because no minima are 
present in the random potential at the microscopic scale 
A. It is important to stress that this random potential is 
superimposed on a mean dimensionless potential u*(</j) 
which itself displays two equivalent minima at the fixed 
point. These global minima, located at the dimension- 
less fields ±<p m ^ (see above), characterize the incipient 
ferromagnetic ordering. Contrary to the disordered elas- 
tic systems for which the random potential dominates at 
large scales over the elastic energy, the random poten- 
tial in random field systems do controls the (sample-to- 
sample) fluctuations, but never becomes predominant in 
the thermodynamic (mean) behavior. The long-distance 
behavior is determined by a subtle interplay between fer- 
romagnetic ordering and randomness. 

Finally, we mention that the rather elusive nature of 
the nonanalyticity in the effective average action may 
become more transparent when studied within a dy- 
namic formalism. As pointed out by Chauve et al.j 4 ^* 
the presence of a cusp may then be related to the exis- 
tence of a nonzero threshold force below which the sys- 
tem stays trapped in a minimum of the random poten- 
tial ("metastable state"). Again, this picture applies a 
zero temperature and it describes an out-of-equilibrium 
dynamic transition, the depinning of an elastic system 
in a random environment^ One may conjecture that 
a similar relation exists in the RFIM between cusp-like 
behavior and zero-temperature driven dynamics among 
metastable states: the analog of the depinning transition 
would then be an "avalanche transition" observed when 
driving the system by slowly varying an applied mag- 
netic field (to sue the terminology of magnetic materials 
in which the phenomenon is commonly observed).— 



B. Role of temperature, droplet phenomenology 
and activated dynamic scaling 

Temperature, as has been stressed in several occa- 
sions in this article and in paper I, plays a peculiar 
role in random field systems. It is irrelevant near the 
fixed points controlling criticality, ordering and quasi- 
ordering. However, it is "dangerously irrelevant" in that, 
at nonzero temperature and slightly off the critical point, 
both static and dynamic quantities display somewhat 
anomalous scaling behavior coming from the scale de- 
pendence on the renormalized temperature! 18 ' 48 To make 
contact with the discussion of the preceding subsection, 
one may summarize the situation as follows ! 42 ' 43 the zero- 
temperature analysis provides information on the typical 
behavior of the system (including the typical fluctuations 
and correlation functions), as described by its ground 
state, whereas small nonzero temperature requires an 
account of rare events such as low-energy excitations. 
Metastable states play a role in both cases: rather in- 
tricate as far as the typical behavior is concerned (see 
above), more direct in the case of low-energy excitations. 



In the latter case, an efficient phenomenological approach 
has been proposed, known as the "droplet picture" i 49 ' 50 
In a nutshell, the droplet approach assumes the exis- 
tence of rare samples (or rare regions in a sample) for 
which, on top of the ground state, an additional mini- 
mum (metastable state) is thermally accessible, having 
an energy above the ground state of the order of the 
temperature. If one defines for a system of linear size L 
sample-dependent (i.e., random) "connected", x c [^], & n d 
"disconnected", susceptibilities as 



Uh] = L- 



(x(x)x(y)) - (x(x))(x(y)) 



and 



Xd[h] = L~ 



ix(x))(x{y)), 



, (41) 



(42) 



where x{ x ) is the fundamental field in the bare action 
and the brackets denote a thermal average for a given 
configuration h of the bare random field (see section II-A 
of paper I). At criticality and at a temperature T (criti- 
cality is attained by fine tuning the bare disorder strength 
whose critical value depends on T as illustrated in Figure 
1 of paper I), most samples, which are characterized by 
a single populated minimum, are such that Xd[h] ~ L 4- * 5 
whereas, due to cancellation of the leading terms in Eq. 
(|4"Tj) . Xc[h] ~ TL 2 ~ V . On the other hand, rare samples, 
which occur with a propability assumed to be of the or- 
der of TL~ e with 8 the temperature exponent given by 
6 = 2 + rj — fj, have Xd[h] ~ Xc[h] ~ L 4 ~'' (since the 
leading contributions to the two terms in Xe[h] no longer 
cancel when two minima are populated). One therefore 
finds for the disorder averaged pth moments: 



X 



P ^ rpjJ){A-fi)-6 



(43) 



(44) 



As a consequence, the fluctuations of the connected sus- 
ceptibility are "anomalous" with, e.g., Xc ^ (Xc) • 

Another important prediction of the droplet approach 
concerns the dynamics of the RFIM near the critical 
point: the critical slowing down of the relaxation is shown 
to be "activated" i 17 ' 18 ' 51 The typical relaxation time di- 
verges exponentially as one approaches the critical point, 
with the effective activation barrier for relaxation diverg- 
ing with system size as 1 ^ L e at criticality.— 

A major step toward formulating a fully consistent 
field theory and renormalization group framework for 
the droplet picture has been accomplished by Balents 
and Le Doussal ) 42 ' 43 ' 45 ' 52 in the context of the random 
elastic model. The core of the connection between FRG 
formalism and droplet phenomenology is the existence 
of a "thermal boundary layer" that governs the highly 
nonuniform limit of the renormalized temperature going 
to zero. In the following, we do not attempt to provide 
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an exhaustive description of thermal boundary layer and 
droplet picture in the RFIM, but rather stress some al- 
ready illustrative results that we obtain within the min- 
imal NP-FRG truncation. 

To study the role of temperature in the NP-FRG for- 
malism, we consider the flow equations for the RFIM, 
keeping now the terms depending on the renormalized 
temperature. From the results of section IV-B of paper 
I, one derives 

- d t u' k (x) = /%r°(x) + 2v d T k l[ d \u' k '(x))u' k »(x), (45) 



-8t8 k (x, y) = Pj=°(x,y) + v d T k { l -[l {d \u'[) + l[ d \u'!_)] 



[l[ d \u'D - ^(u'L)]S^\x,y)] + + tfV-)] 



M/„'MU( 01 )/ 



x C\x,y)] - m a >(u';)+l?>(u'L)]Sr(x,y) 



(46) 



where (3^i =0 and /3j =0 are the T = beta functionals 
given by the right-hand sides of Eqs. (f22| and (JTTJ) re- 
spectively in which the running exponents r\ k and fj k are 
now expressed at T ^ 0, and where u± = u k (x ± y) as in 
Eq. (TTTl). 

Taking the limit y — > in Eq. (|4"6")l and allowing for 
cusp-like behavior, one finds for 5 k fi{x) — S k (x 1 y = 0) 
that 

-d t S kfi (x) = p So \ reg {x) 

- y4 d) «W) d 2 y (6 k (x,y) - 4,o(*)) 2 | y=0 (47) 
+ v d T k l[ d \u' k , (x))d 2 6 k (x,y)\ v=0 

with 

Ps \reg(x) = Pj = °\reg{x) + V d T k 



l i d \u'k(x))S k0 (x) 



2lf\u'l{x))5' kfi {x)u'l!{x) 



(48) 



and f3j~°\reg is given by the right-hand side of Eq. (fT9l . 

If a cusp appears in the renormalized cumulant S k (x, y) 
near the fixed point, i.e., 5 k (x,y) = 5 k ,o(x) + \y\S k , a (x) + 
0(y 2 ) as in Eq. then the two "anomalous" terms 

appearing in the right-hand side of Eq. (|47|) behave quite 
differently: the term characteristic of the T — behav- 
ior, d 2 (5 k (x,y)-5 kfi {x)) 2 \ v=0 , goes to 25» ja (x) 2 , whereas 
that proportional to T k , d 2 5 k (x, y)\y=o, blows up. For a 
fixed point to be reached and the theory be renormaliz- 
ablc, the latter divergence must be cancelled. The solu- 
tion is that (i) there should be no cusp at finite T k and 
(ii) convergence to the cuspy T = fixed-point function 
is nonuniform in y as T k — > and takes the form of a 
boundary layer. 



In the close vicinity of the fixed point, when both y and 
T k approach zero, one anticipates the following behavior: 

5 k (x,y) = S., (x) + T k f(x, |-) + 0(T 2 ) (49) 

where 5*.o(x) is the (T = 0) fixed point result for y = 
and f(x, y) is a scaling function which is even in x and y 
and analytic around y = 0; 0(T k ) denotes terms of order 
at least T k at fixed scaling variable y. With Eq. (|49|) . 
the two anomalous contributions described above can be 
expressed as 

d 2 y (6 k (x,y)-5 kfi (x)) 2 \ y=0 = 2f(x,0)f(° 2 \x,0), (50) 

T k d 2 6 k (x,y)\ y= o = f(° 2 Hx,0), (51) 

with the same convention as before for the notation of 
partial derivatives. 

Inserting the above expressions into the flow equation 
for 5 k (x,y), Eq. ([46]) . and using the fact that S kt o(x) is 
solution of the the fixed-point equation at zero tempera- 
ture for y = 0, one finds up to a 0(T k ) that the scaling 
function / must satisfy that 



I/CO/,..", .', , :>~! 

2 2 



(u'tix)) d*(f(x,y) - f(x,Q)Y 



(52) 



is independent of y. 

The solution is easily obtained as 

f(x,y)-f(x,0) = 
l { x\<{x)) 



l {d \K{x)) 



\ 



' lt\<{x))f^){x,Q)) 
l { '\<(x)) 



!J 2 



(53) 

with f( 02 >(x,0) < and u"(x) given by the fixed-point 
solution of Eq. (|4"5")) . The zero-temperature cuspy fixed 
point is recovered by considering y ^ and T k — > 0, 
which leads to 



f(n,y-> ±°°) ~ 5*,a( x )\y\ 



(54) 



with <5*, (x) = /( 02 )(x,0) < 0. On the other hand, when 
y — > at fixed T k ^ 0, i.e., near y = 0, f(x, y) — f(x, 0) = 
0{y 2 ). Eq. (f53|) thus describes the rounding of the cusp 
near y = in a layer whose width of order T k goes to 
zero as k — * 0. 

We can now make contact with the droplet description 
of the RFIM. From the effective average action T k , one 
has access to all Green functions of the system via the 
1 — PI vertices (see paper I). In the present minimal 
truncation, 

[{0a}] = / {^E [\Zm,k\d4> a {*)\ 2 + Wa(*))] 
-^2 E V k ((t) a (x),(t) b {x)) 



a, 6=1 



(55) 
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where we use an explicit dependence on a bare tempera- 
ture T for book-keeping purpose. 

The 1 — PI vertices are obtained by functional dif- 
ferentiation and their expression is given in Appendix 
B. When all fields are taken as uniform and, moreover, 
equal, the first vertices have the following form: 

{<y„&f< 2) (0 ;gi )+ff (<M;gi)}, (56) 



than the running scale k. Thanks to the RG framework, 
this is enough to provide a determination of the anoma- 
lous dimension of the field 77, but in what follows, we only 
consider the case of external momenta set to zero. 

Next, focusing on the critical (scaling) region, we in- 
troduce dimensionless functions and fields by using the 
scaling dimensions suitable for a zero-temperature fixed 
point (see section II); so, for instance, 

~E#(0) ~ (62) 



r 



(3) 

k,(a,q 1 )(b,q 2 )(c,q 3 )({ < t ) f = </>}) 

(2n) d S( qi 



92 



-93)lfL(0), 



r (4) fe,(Q,q 1 )(b,q 2 )( c ,q 3 )(d,q 4 )({ ( ? !) / = 4>}) = 

(2n) d S( qi + g 2 + 93 + q^ abc M, 



(57) 



(58) 



where the functions ff (cf>;q), faftq), T^tf), 

T ( ^ abcd (4>) can be derived from Eqs. (|BlHB3|) . 

The (connected) Green functions may be obtained 
from the replicated free energy functional Wfe[{J }], 
which is the Legendre transform of the effective aver- 
age action r fc [{^ a }]. The procedure consists of using the 

standard formulas that relate the Wjf''s to the T^'s^ 
and expanding both sides in number of free replica sums 
as explained in paper I, keeping only the leading terms. 
This is detailed in Appendix B. The Green functions can 
be cast in a form similar to that of the 1 — PI vertices, 
namely, for equal field arguments, 



(2), 



,(3) 



w 



k,(a,qi)(b.q2) 
(2) 



4>}) = (2n) d 5( qi + q 2 )x 

(2) M ^0 (™> 



S ab G\f> (0; qi ) + (0, <j>;qx) + 0(]T) 

/ 



W, 



(3) 



fc,(Q,qi)(6,q 2 )( C ,q 3 )(^/ ~ 



92 + g 3 )G^ bc (0; 91,92), 
(60) 



W 



(4) 

k,(a,qi)(b,q 2 )(c,q 3 )(d,q 4 ) 



({4>f = 4>}) 

(4) 



(2-K) a S{ qx + g 2 + g 3 + qi)Gl ' abcd ((t>; gi,g 2 ,93), 



(61) 

where the various quantities appearing in the right-hand 
sides are related to the 1 — PI counterparts in Eqs. ((561 
I5"5f as discussed in Appendix B. Note that in the present 
truncation of the effective average action (limited to the 
first order of the derivative expansion), 13 information on 
the Green functions is essentially limited to zero external 
momenta, more precisely external momenta with I9I less 



fc 4 ~ rl 

-=2-&(<Pi,¥>2)- (63) 



In addition, we eventually take the limit of equal field 
arguments in all expressions after insertion of the results 
obtained for the thermal boundary layer description of 
5 k (<pi,<p2), Eqs. (1491531) . 

In the scaling region, the Green functions at zero ex- 
ternal momenta can then be expressed as 



Gi 2) (</»,0)~T fe fc-( 4 -^/il 2) (^), 



(2), 



,«ko)~Ar( 4 -SrV), 



(64) 



(65) 



G 



(3) 



0,0) ~ k-i-^-V\gi 3 \<p) + T k h? ) (tp)(6 ab 



+ S bc + S ca ) + 0(Ti) 



(66) 



G 



(4) / 
k.abcd \ 



; 0,0,0) ~ k- d - 2 ^\g^^) + T k h?\cp)x 



(S ab + S ac + 5 a d + S bc + S bd + 5 cd ) + t J fl ^^ 



0(Ti) 



(67) 



where — > (as k^ d i+ ^/ 2 ^>Lp) and where the fonctions 

gi p \ip), hi p \ip), p — 2,3,4, ... are obtained from u"(tp), 
5* t o(ip) and their derivatives. Their expression is not par- 
ticularly illuminating and we do not reproduce them here; 
see Appendix B for more details. 

From the definition of the replicated generating func- 
tional Wfc[{J a }], one can derive the relation between the 
replica Green functions considered above and the phys- 
ical Green functions directly defined in the disordered 
system. 13 Still working at the running scale k and at 
zero external momenta, the first moments of the (ran- 
dom) "connected" and "disconnected" susceptibilities in- 
troduced in Eqs. (|41l42p are for instance given by 



Xk.c^G^^-,0) 



(68) 
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xm = 4 2 W; o), 

whereas the second moments read 

,(4) 



(4) 

k.abcd 



Xk,d 



k d G<Cl bc Mi 



(69) 



(*)], (70) 



(71) 



where distinct replica indices here mean distinct replicas 
(no summation implied). 

Putting together all the above results, Eqs. (|64H7ip . 
we find that the moments of the random "disconnected" 
suceptibility scale as 



2 ~ 1.-2(4-9) 



Xfc,d 



(72) 



(73) 



whereas those of the random "connected" susceptibility 
scale as 

JTd-T k k-^ ~Tk-^\ (74) 



r ,- 2 (4-^ / (02) (^0) 

<W (75) 

y fc -2(4-rl)+e 



where we have used that Xfc ~ T/c e , and we recall that 
f {02) (f, 0) = £" a (<p) < 0. Notice that only the term due 
to the boundary layer appears in Eq. (|75p (the other 
contributions cancel out). 

From this analysis, we therefore obtain that the mo- 
ments of the random "connected" and "disconnected" 
susceptibilities in the truncated NP-FRG precisely scale 
as in the droplet description: compare Eqs. (|72ll75p and 
Eqs. (|43l44p . with L — fc _1 . The "anomalous" scaling of 
the moments of the "connected" susceptibility, which is 
due to rare low-energy excitations in the droplet picture, 
results in the NP-FRG from the presence of a thermal 
boundary layer in the vicinity of the zero-temperature 
fixed point, as illustrated by Eq. ([75"|). This is in complete 
agreement with the more detailed analysis performed in 
Ref. [HI for disordered elastic systems. 

To conclude this section, we briefly address the ques- 
tion of the slowing down of the relaxation toward equilib- 
rium near the critical point. At long times, the dynamics 
of the RFIM can be modeled by a Langevin equation, 



d T x{x, 



-n 



SS[ X ;h] 



(76) 



where r denotes the physical time (to be distinguished 
from the RG "time" t) and ((x, r) is a thermal noise 
taken with a gaussian distribution characterized by a zero 
mean and a second moment 



In Eq. ([76]) S[x; h] is the bare action for the RFIM, with h 
being the bare random field (see for instance section II of 
paper I) , and fHs a kinetic coefficient that describes the 
bare relaxation rate and sets the elementary time scale 
in the problem. (We consider here the case of a non- 
conserved order parameter^ but the case of a conserved 
order parameter could also be of interest^) 

An RG formalism can be conveniently implemented 
by using standard field theoretical techniques to build 
the generating functional of the time-dependent corre- 
lation and response functionsi 53 i 56 ' 57 i 58 Associated with 
this functional is a "bare dynamic action" that depends 
on two fields, the fundamental field and a "response" 
field. The average over the quenched disorder can now 
be performed without introducing replicas. By a Legen- 
dre transform, one then introduces an "effective dynamic 
action" which is the generating functional of the disorder- 
averaged, time-dependent 1 — PI vertices. 

In this setting, one can repeat the steps detailed in pa- 
per I to construct a NP-FRG approach to the dynamics: 
add a mass-like regulator with an infrared cutoff function 
that suppresses the contribution of low-momentum and 
low frequency modes, define a (dynamic) effective average 
action at scale k, whose evolution with k is governed by 
an exact RG flow equation, devise a truncation schemed 
In the present problem, this latter step can be done in 
the spirit of the minimal trucation considered above and 
in paper I. One however needs an additional assump- 
tion concerning the time dependence which, similarly to 
the spatial dependence, can be handled via an appropri- 
ate "derivative expansion" . The simplest approximation 
that captures the physics is a "single time scale" approx- 
imation in which one introduces a single renormalized re- 
laxation rate Q k - (This parallels the single wavefunction 
renormalization parameter used to describe the spatial 
dependence of the field in the minimal truncation: see 
paper I.) 

From the boundary layer structure and by analogy 
with the previous work on the random elastic model ; 41 ! 43 
one then expects that the renormalized relaxation rate 
flows as 



d t ln(ft k ) 



-t: 



(78) 



C(as, r)C(y, r') = 2T05(r - t')6(x - y) 



(77) 



near the zero-temperature fixed point. This indeed cor- 
responds to activated dynamic scaling and fits in with 
the droplet picture summarized above. A proper deriva- 
tion of this result and an account of the (expected) broad 
distribution of relaxation rates would require a detailed 
dynamic treatment, but this goes beyond the scope of 
the present article. 



VI. CONCLUSION AND PERSPECTIVES 

In this paper, which is the second part of a series of ar- 
ticles reporting our work on a nonperturbative functional 
renormalization group (NP-FRG) approach for random 
field models and related disordered systems, we have 
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applied the formalism presented in paper I to the d- 
dimensional random field O(N) model. We have focused 
on two main issues related to the long-distance physics of 
the model: the breakdown of the dimensional reduction 
property predicted by conventional perturbation theory 
and the nature of the phase diagram and ordering tran- 
sitions in the (N, d) plane. 

Within our NP-FRG approach, the way out of dimen- 
sional reduction is the appearance of a strong enough 
nonanalyticity in the field dependence of the dimen- 
sionless effective average action near the relevant zero- 
temperature fixed point. We have shown that this oc- 
curs below a critical dimension dj)n(N), which goes con- 
tinuously from N DR = 18 as d — ► 4+ to d^R — 5 for 
N = 1. In addition, we provide a description of critical- 
ity, ferromagnetic ordering, and quasi-long range order 
in the whole (N, d) plane. The NP-FRG method is able 
to directly address the phase diagram of the model in 
low (physical) dimension d and small (physical) number 
of components N: in particular, we find that there is no 
"Bragg-glass" phase, i.e., no phase with quasi-long range 
order in the 3-dimensional RFAFM. Note that all those 
results are made possible by the very structure of the 
present RG approach which is (1) functional, (2) approx- 
imate but nonperturbative, and (3) devised to provide a 
continuous and consistent description of the whole plane 
of N, d. 

Building upon earlier work on random elastic 
modelS ) 42 i 43 we have also shown how the NP-FRG for- 
malism gives access to both the typical behavior of the 
system, controlled by zero-temperature fixed points with 
a nonanalytic dimensionless effective action, and to the 
physics of rare low-energy excitations ("droplets"), de- 
scribed at nonzero temperature by the rounding of the 
nonanalyticity in a thermal boundary layer. 

Work still remains to be done for a complete under- 
standing of random field models. We have pointed in 
several occasions in this paper and in the preceding one 
that clarifying the putative link between breaking of the 
underlying supersymmetry^ and appearance of a nonana- 
lyticity in the dimensionless effective action would require 
to "upgrade" the present NP-FRG formalism to a super- 
field formulation of the random field models. We defer 
this, as well as the study of improved nonperturbative 
truncations, to a forthcoming publication. We have also 
indicated an interesting extension of the present work to 
the dynamics of the random field Ising model (RFIM), 
both to the out-of equilibrium driven dynamics at zero 
temperature and to the (activated) relaxation to equi- 
librium at nonzero temperature. Finally, the connection 
to the proposed picture of the RFIM in terms of "spon- 
taneous replica symmetry breaking" and "replica bound 
states" 60 ' 61 i 62 ' 63 remain to be investigated. 

The NP-FRG formalism appears as a powerful tool to 
study random field models and related disordered sys- 
tems. Whether such an approach can be generalized 
to tackle another major unsettled problems of the field 
of disordered systems, the long-distance physics of spin 



glasses, is a challenging but completely open question. 
We thank D. Mouhanna for helpful discussions. 



APPENDIX A: NONANALYTICITY IN THE 
RFIM NEAR d = 6 

Our starting point is the RG flow equations for 6k,o(x) 
and 5k,2p(%) obtained by assuming that Sk(x,y) is reg- 
ular enough near y = 0: see Eqs. (|19II21[) . The linear 
operators appearing in those equations are given by 

L<2p[u", So, 8 2 ] = - [p(d -4 + fj k ) + 2r] k - fj k ] - 

i(d - 4 + fj k )xd x + 2vJl {d \u"(x))5 (x)d 2 +2(p+ 1) 



l 2 d \u"(x))S' (x) - 2li d \u"( x ))5 ( x )u'l'(x) 



p(2p + 3) + 1 



li d \u"(xM(x)-2li d) (u"(x))5 (x)x 



u'£(x) + 2li d \u"(x))5o(x)u'>'(x) 2 



l {d \u"{x))52{x) 



P{2p + i). 



(Al) 



for p > 2; the expression of L2[u",Sq] is obtained by 
setting p — 1 in the above equation and dropping the 
last term so that S2(x) no longer appears in the operator. 
Note that, as for the threshold functions (see paper I), 
there is an explicit dependence on k due to r]k and fjk that 
comes on top of the dependence that may occur through 
the arguments. 

Near d = 6, one finds, as developed in section V- 
A of paper I. that the fixed point is characterized by 

= f h = 0(e 2 ), u»{x) = e(\^/2)(3x 2 - x^) + 0(e 2 ), 

S*(x,y) = 1 + e 2 d(x,y), with x 2 m ^ = 6v e l ( 2 \o), Ai* = 

(36v 6 4 6) (0)r\ and d{x,y) = 0(1). The result for 
5*(x,y) implies that (x) = 0(e 2 ) and 5* 2 P (x) = 
0(e 2 ). 

After inserting these results in Eq. (|A1[) . one obtains 
that for p = 0(1), 

L 2p * [x) ~~2p-xd x + 2v d l (d) (0)d 2 + 0(e) (A2) 

at the fixed point. We have made the implicit assump- 
tion, whose consistency can be checked, that the deriva- 
tives with respect to x acting on the 5k,2p's do not modify 
the order in e. The eigcnfunctions of the above linear op- 
erator (with the condition that they are bounded by poly- 
nomials at large values of the argument^! are the Hermite 

polynomials H n (x/\J Avdl 2 d \^)^r^ with associated eigen- 
values \2 P ,n — —(2p + n), with n 6 IN. Recalling that the 
RG time t in Eq. (f2T|) goes to —00 as k — * 0, the above 
result means that the corresponding directions are irrel- 
evant on approaching to the fixed point. (This confirms 
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the result found in section V-A of paper I that the fixed 
point given above is once unstable at first order in e.) 

However, a new phenomenon may appear when p is 
very large and scales as 1/e 2 . In this case, one finds that 

L 2p *(x) ~ -2p jl - V -ip[2lt\Q)u'i'{xf + 4 d) (0)x 

(<5^)-Mx))] + 0(e)j, 
(A3) 

where the whole second term in the braces is of order 1 
and, if positive, can become larger than 1 for some value 
of p so that L2p* becomes positive. 

To analyze the sign of I-2 P *, one has to study Sq^(x) 
and 82* (x) at order e 2 . This is easily performed from 
Eqs. ([15 ]) and ([20 ]) . One finds that S'^(x) = -6 2 *(x) = 
lSi^^WA^e 2 , which indeed guarantees that the sec- 
ond term in the right-hand side of Eq. (|A3|) is positive. 
Evaluating L 2p * for x = x m * and using the expression of 
u'"(x) now gives 

L 2p *(x 7W ) ~ -2 Jl - K 2 (pe 2 ) + 0(e)), (A4) 



with K = 4 (0)4 (0) /(6Z3 (0) ) 2 . One therefore con- 
cludes that forp > l/(i^e) 2 L2 P *{x m *) becomes positive, 
which, according to Eq. (f2Tj) , leads to a divergence of 
5k,2p{xm*) as k — ► 0. As a consequence, the renormal- 
ized disorder cumulant displays a subcusp of order 1/e 2 
at the fixed point. A related phenomenon has also been 
proposed by Feldman. 25 



{^C(^) - ^[S abc A {20 \<p c ,cb d ) +perm(abcd)} 

- 7^ [SabScdA^ {4>a, 4>c) + SacSbdA^ (4>a, <j>d) 

(B3) 

where S a bc = 5 ab &bc, S a bcd = S a bS bc S c d, a perm(abc)" de- 
notes the two terms obtained by circular permutations 
of the indices a&c, and "perm(a6cd)" denotes the three 
terms obtained by circular permutations of the indices 
abed. All other notations are as in paper I and above. 

The connected Green functions 

Wj?) , , >({<Pf}) are related to the 1 - PI 
vertices by formulas deriving from the Legendre trans- 
form between W k and For instance, the 2-point 

(2) 

connected Green function Wj. is the inverse of the 

2-point 1 - PI vertex, = rf )_1 . By using the 

expansion in number of free replica sums detailed in 
section II-D of paper I, one obtains at leading order 

<t, qi )(6, g2 )({M) =(^) d S(Qi + <Z2){<UGf (4> a ; qi ) 

+ Gi 2) (0 a ,0 b ;< Zl ) + O(^)], 

(B4) 

where O(J^) denotes higher orders in the expansion in 
number of free replica sums and with 

a " to) ° wfw > (B5 > 



APPENDIX B: GREEN FUNCTIONS IN THE 
TRUNCATED NP-FRG OF THE RFIM 



The 1 — PI vertices are obtained by functional differ- 
entiation of the truncated effective average action given 
in Eq. (|55p . For uniform field configurations one finds 



.(2) 



y ab ^;[Z m>k ql + U' k ' {(j)a)] - ^Afe(0 a ,, 



(Bl) 



■ perm(abc)] 

(B2) 



G 



(2)/ 



<>;q) = 



A fc (0 a ,0fc) 

Z m . k q 2 + UH{<f> a ))(Z m<k q2 + ' 

(B6) 

The 3- and 4-point connected Green functions are de- 
rived along similar lines, using the standard graphical 
representation^ and keeping the lowest order in the ex- 
pansion in free replica sums. 

Consider now the scaling region. As discussed in the 
main text around Eqs. ([521 [53J), one can introduce di- 
mensionless quantities and use the results concerning the 
thermal boundary layer. For the I— PI vertices evaluated 
at zero external momenta and for equal field arguments, 
one explicitly obtains the expressions of the functions 



=i(2) p(2) r (3) 



k,abc ' 



^(2) 



k,abcd 



appearing in Eqs. (|56ti58|) : 



kq = 0) 



(B7) 



ri 2) (0,0;g = O) 



.^-S,, (<p), (B8) 
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V k'abM) ~ f <<W** {<P)- 

k Inserting the above equations in the expressions of the 

^*,o(W ^ g si connected Green functions, Eqs. (|B4BB6|1 and their gen- 

2Xfc bc ca J ' eralizations for higher-order functions, finally leads to 

(B9) Eqs. (|64ti67[) of the text, with for instance 

(4) k 4-n-(d-4+v) j fm ^ *i' j0 (^) rj ffi 2) M = *.,oM<'M" 2 (Bll) 



r fc,a&cd(^) - jT |<Ucd<"(<?) - * 4Tfc [<Uc 

+ Sbcd + Soda + Sdab + Sabred + S ac Sbd + SadSbc] 

/ (02) (V,0) r, ^ h<?\<p)=u , :(<p)- 1 > (B12) 

7777? "abc + Obcd + O c da + Odab 

4T 2 L 

., 1 etc... As stated in the text, the expressions for the other 

— (Oa&Ocd + OacO&d + OadOfcc) f , r (p) / \ u {p) , \ n a 1 • 1 II f \ 

J J functions^; (<p), hi (<p), p = 3, 4,..., only involve u"(ip), 

(BIO) <5*,o(^) and their derivatives. Obtaining them is tedious 
but straightforward, and the resulting formulas are not 
where of course the dimensionful field (f> vanishes (as worth displaying. 
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